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Abstract 


This  report  documents  the  usage  of  the  routine  FITS,  which  provides  au¬ 
tomated  fits  of  various  analytical,  commonly  used  probability  models  from 
input  data. 

This  routine  is  intended  to  complement  the  previously  distributed  rou¬ 
tine,  FITTING,  documented  in  RMS  Report  14  (Winterstein  et  al,  1994). 
The  FITTING  routine  implements  relatively  complex,  four-moment  distribu¬ 
tion  models,  whose  parameters  are  fit  with  numerical  optimization  routines. 
While  these  four-moment  fits  can  be  quite  useful  and  faithful  to  the  observed 
data,  their  complexity  can  make  them  difficult  to  automate  within  standard 
fitting  algorithms. 

In  contrast,  the  routine  FITS  is  intended  to  provide  more  robust  (lower 
moment)  fits  of  simpler,  more  conventional  distribution  forms.  For  each 
database  of  interest,  the  routine  estimates  the  distribution  of  annual  max¬ 
imum  response,  based  on  the  data  values  and  the  duration,  T,  over  which 
they  were  recorded.  To  focus  on  upper  tails  of  interest,  the  user  can  also 
supply  an  arbitrary  lower-bound  threshold,  above  which  a  shifted  dis¬ 
tribution  model — exponential  or  Weibull — is  fit.  (In  estimating  the  annual 
maximum  response,  the  program  automatically  adjusts  for  the  decreasing 
rate  of  response  events  as  the  threshold  xiow  is  raised.) 

Note  also  that  FITS  permits  modelling  of  multiple  response  databases 
which  may  represent  different  statistical  populations.  Examples  include  (1) 
wave  or  response  statistics  due  to  waves  from  different  directions,  or  (2)  re¬ 
sponses  of  Gulf  of  Mexico  structures  due  to  hurricanes  alone,  eddy  currents 
alone,  or  combined  hurricane-eddy  events.  By  including  the  relative  occur¬ 
rence  rates  and  response  distributions  for  each  population,  FITS  estimates 
the  total  contribution  to  the  annual  maximum  response. 
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1  Distribution  Fitting:  General  Methodol¬ 
ogy 


This  report  documents  the  usage  of  the  routine  FITS,  which  provides  auto¬ 
mated  fits  of  various  analytical,  commonly  used  probability  models  from  in¬ 
put  data.  This  routine  is  intended  to  complement  the  previously  distributed 
routine,  FITTING,  documented  in  RMS  Report  14  (Winterstein  et  al,  1994). 
The  FITTING  routine  implements  relatively  complex,  four-moment  distribu¬ 
tion  models,  whose  parameters  are  fit  with  numerical  optimization  routines. 
While  these  four-moment  fits  can  be  quite  useful  and  faithful  to  the  observed 
data,  their  complexity  can  make  them  difficult  to  automate  within  standard 
fitting  algorithms. 

In  contrast,  the  routine  FITS  is  intended  to  provide  more  robust  (lower 
moment)  fits  of  simpler,  more  conventional  distribution  forms.  For  each 
database  of  interest,  the  routine  estimates  the  distribution  of  annual  max¬ 
imum  response,  based  on  the  data  values  and  the  duration,  T,  over  which 
they  were  recorded.  It  can  also  estimate  the  annual  maximum  response 
distribution  due  to  multiple  types  of  response  “events,”  intended  to  reflect 
different  statistical  populations.  Examples  include  waves  from  different  di¬ 
rections,  or  response  events  due  to  hurricanes,  eddy  currents,  or  combined 
hurricane-current  events. 

Specifically,  FITS  first  estimates  the  individual  response  CDFs  Fi(x)=P[ 
Outcome  <  a:  in  database  i]  from  the  data  in  each  database.  It  then  predicts 
the  corresponding  distribution  of  annual  maximum  response  by  assuming 
that  events  from  each  database  occur  in  independent,  Poisson  fashion: 

P[  Annual  max  >  x]  =  Gann{x)  =  1  -  exp[-  ^i/,C?t(x)] ;  G,(x)  =  1  —  Fi(x) 

(1) 

Here  Vi  is  the  mean  annual  rate  of  events  in  database  z,  estimated  from  the 
observed  number  per  year  (possibly  over  a  user-defined  threshold;  see  below). 
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1.1  Available  Distribution  Types 


Specific  distributions  currently  included  in  FITS  to  estimate  Fi{x)  include 
the  following,  as  catalogued  by  the  distribution  index  IDIST; 


IDIST=1: 

Normal  Distribution 

IDIST=2: 

Lognormal  Distribution 

IDIST=3: 

Exponential  Distribution 

IDIST=4: 

Weibull  Distribution 

IDIST=5: 

Gumbel  Distribution 

IDIST=6: 

Shifted  Exponential  Distribution 

IDIST=:7: 

Shifted  Weibull  Distribution 

The  distributions  IDIST=1  through  5  are  all  fit  to  statistical  moments  of  all 
available  data.  The  single-parameter  exponential  preserves  only  the  mean 
nix  of  the  data,  while  the  normal,  lognormal,  Weibull,  and  Gumbel  preserve 
both  the  mean  and  standard  deviation  estimated  from  the  data.  The 
shifted  versions  of  the  exponential  and  Weibull  distributions  (IDIST=6  and 
7)  employ  a  user-defined  lower  threshold  xiow,  and  fit  standard  exponen¬ 
tial/ Weibull  models  to  rr  -  based  on  observed  moments.  (In  estimating 
the  annual  maximum  response,  the  program  automatically  adjusts  for  the 
decreasing  rate  of  response  events  as  the  threshold  xiow  is  raised.) 

The  result  aims  to  provide  the  user  with  a  suite  of  smooth  probability 
models,  to  be  fit  throughout  the  body  of  the  available  data.  It  does  not 
directly  address  various  special  topics  of  data  fitting;  e.g.,  selective  tail  fit¬ 
ting,  fitting  bimodal  models  to  hybrid  data,  etc.  Some  of  these  issues  can 
be  addressed,  in  a  limited  way,  through  the  use  here  of  the  shifted  models 
(IDIST=6  and  7).  In  this  way  the  user  can  focus  the  distribution  modelling 
resources  on  the  extreme  response  levels  of  interest. 
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More  specific  tail-fitting  procedures  have  not  been  given  here,  because 
optimal  use  of  these  tends  may  be  rather  problem-specific.  In  the  same  vein 
our  extremal  models  are  limited  here  to  so-called  “Type  F  behavior,  leading 
to  (shifted)  exponential  distributions  of  peaks  over  a  given  threshold  and 
to  Gumbel  distributions  of  annual  maxima.  Type  II  and  III  distributions 
are  ill-suited  to  our  moment-fits,  due  to  potential  moment  divergence  (Type 
II)  or  to  the  difficulty  in  predicting  truncated  distributions  (Type  III)  from 
moment  information. 


1.2  Limitations 


Several  parameters  have  been  assigned  maximum  values  in  the  routine  FITS. 
These  include  the  upper  limits 


•  NFMAX,  the  maximum  number  of  files  (databases)  to  be  fit,  has  been  set 
to  3. 

•  NMAX,  the  maximum  number  of  data  per  database,  has  been  set  to 
32768. 


Both  of  these  limits  have  been  set  in  PARAMETER  statements  in  the  main 
driver  program  to  FITS.  These  are  rather  arbitrarily  selected  limits,  and  can 
be  reset  by  the  user  without  fundamental  consequence. 
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2  Distribution  Fitting:  Routines 


The  fitting  algorithm  calls  the  following  set  of  subroutines: 


CALMOM:  Estimates  the  mean  rrix,  standard  deviation  Ox,  skewness  03  and 
kurtosis  04  from  an  input  set  of  data.  These  are  based  on  unbiased 
estimates  of  the  cumulants  ki=mx,  ^2=0^,  k3=a3<rl,  and  fc4=(a4— 

If  the  user  includes  an  optional  lower  limit  Xminj  moments  of  the  shifted 
variable  (a:  -  x,„in)'^=max(0,a:  -  Xmm)  are  estimated. 

DISPAR:  Based  on  the  sample  moments  estimated  in  CALMOM,  DISPAR  seeks 
a  consistent  set  of  distribution  parameters.  The  interpretation  of  these 
parameters  depends  on  the  distribution  type  selected  by  the  user.  Ap¬ 
pendix  A  includes  a  complete  listing  of  the  distribution  functions  and 
their  parameters. 

GETCDF;  For  the  user-defined  distribution  type  with  the  distribution  param¬ 
eters  from  DISPAR,  this  routines  estimates  the  cumulative  distribution 
function  value,  F(a;)=P[Outcome  <  x]  for  given  input  x  value. 

FRACTL:  For  the  user-defined  distribution  type  with  the  distribution  param¬ 
eters  from  DISPAR,  this  routines  estimates  the  fractile  x  corresponding 
to  a  specified  input  value  of  the  probability  p=F(a:)=P[Outcome  <  x]. 

QDMOM:  Uses  Gaussian  quadrature  to  estimate  the  first  four  moment  of  the 
theoretical  fitted  distribution.  These  can  be  compared  with  the  sample 
moments  from  the  data,  as  given  by  CALMOM,  to  verify  the  accuracy  of 
the  fitted  model— and  in  the  case  of  the  higher  moments  not  used  in 
the  original  fitting,  to  test  its  accuracy. 


The  routines  GETCDF  and  FRACTL,  which  supply  general  distribution  func¬ 
tions  and  their  inverses,  may  also  be  useful  in  other  stand-alone  applications; 
e.g.,  to  create  a  distribution  library  for  standard  FORM/SORM  or  simulation 
analyses  (Madsen  et  al,  1986),  or  for  use  with  new  Inverse  FORM  algorithms 
(Ude  and  Winterstein,  1996). 
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3  Input  Format  and  Wave  Height  Example 


3.1  Data  Input 


The  file(s)  containing  data  are  read  in  free  format,  one  datum  per  line.  Non- 
numerical  input  are  taken  as  comments  and  ignored.  The  first  numerical 
value  found  is  taken  to  be  the  duration  of  the  database  (in  years).  Remaining 
values  (1  per  line)  are  interpreted  as  data,  and  data  are  read  until  the  end 
of  the  file  is  encountered. 

We  illustrate  the  input  here  through  a  simple  example,  involving  signifi¬ 
cant  wave  height  data.  These  are  in  fact  19  annual  maximum  Hs  values,  esti¬ 
mated  by  hindcast  in  a  Southern  North  Sea  location  (Winterstein  and  Haver, 
1991).  These  data  are  stored  in  appropriate  format  in  the  file  gumbel.dat. 

The  contents  of  gumbel .  dat  are  listed  below.  The  first  line  reflects  that 
19  years  are  covered,  and  the  remaining  19  lines  contain  the  actual  maximum 
Hs  encountered  in  each  of  these  19  years.  (In  this  input  file  the  data  are 
given  in  descending  order;  this  is  not  required  by  the  program.) 


19. 

9.66 

9.44 

9.18 

9.17 

8.85 

8.79 

8.60 

8.58 

8.54 

8.49 

8.09 

8.08 

8.06 

7.47 
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7.42 

7.41 

7.31 

7.16 

6.92 


Contents  of  gumbel.dat. 


3.2  Runtime  Input:  Batch  Mode 


We  seek  to  invoke  FITS  under  the  following  conditions: 


1.  Results  to  be  written  to  a  file  named  gumbell.out.  (The  distinction 
between  lower-  and  upper-case  letters  in  filenames  is  honored  by  Unix, 
and  ignored  in  DOS.) 

2.  Distribution  results  are  to  be  written  for  x  (wave  height)  values  ranging 
from  XMIN.=5.0  to  XMAX.=20.0  m,  at  increments  of  DX.=0.5  m. 

3.  There  is  only  NFILES=1  database  to  be  fit. 

4.  The  wave  height  data  are  stored  in  a  file  named  guinbel.dat. 

5.  The  user  desires  to  fit  a  Gumbel  distribution  (IDIST=5)  to  these  data. 


In  this  example  FITS  requires  5  lines  of  input,  corresponding  to  the  in¬ 
formation  given  in  the  5  items  above.  In  this  case  the  5  input  lines  are  as 
follows: 


gumbell.out 
5.0  20.0  0.5 

1 

gumbel . dat 
5 


Input  line  1: 
Input  line  2: 
Input  line  3: 
Input  line  4: 
Input  line  5: 


File  where  output  is  to  be  written 
xmin,  xmax,  dx  for  writing  dist  results 
number  of  databases  (input  files) 
name  of  1st  database  (input  file) 
idist  =  dist  type  (idist=5  ->  Gumbel) 
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{the  and  following  information  is  not  read  by  FITS;  these  are  given  here 
to  remind  the  user  of  the  input  definition).  These  5  lines  have  been  stored  in 
a  file  named  gumbell .  in.  Thus,  to  execute  FITS  in  a  batch  mode,  the  user 
can  type  (in  either  Unix  or  DOS  environment)  the  following  command: 


fits  <  gumbell. in 


This  will  create  (or  overwrite)  a  corresponding  file  gumbell. out,  whose 
contents  is  discussed  in  the  next  section.  While  FITS  is  executing,  the  user 
will  see  prompts  for  terminal  inputs;  these  can  be  simply  ignored  (or  directed 
toward  null  device)  in  this  batch  mode  operation. 
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3.3  Runtime  Input:  Interactive  Mode 


If  the  user  simply  types  “FITS”,  he  will  be  prompted  for  each  input  (the 
same  5  lines  in  this  case)  with  interactive  explanations.  For  example,  be¬ 
fore  requesting  the  distribution  index  IDIST  the  program  lists  all  available 
distribution  types  and  associated  IDIST  values.  Inputs  with  invalid  formats 
are  ignored.  This  interactive  mode  may  be  particularly  useful  for  first-time 
users.  (The  text  with  input  prompts  is  written  to  the  logical  unit  lOERR, 
which  is  set  to  0  in  the  driver  program  for  FITS.  The  user  can  reset  this  if 
necessary.) 

The  following  is  a  screen  dump  of  the  terminal  input  prompts  and  user’s 
response.  Lines  beginning  with  “>”  are  input  prompts  generated  by  the 
program.  Other  lines  are  the  user’s  response;  in  this  case  there  are  precisely 
5  lines  of  response,  as  in  the  batch  mode  input  file  given  previously. 

> 

>  ♦♦  ENTER  FILENAME  WHERE  OUTPUT  WILL  BE  WRITTEN  ♦♦ 

> 

>  ENTER  OUTPUT  FILENAME: 

gumbell.out 

> 

>  *♦  ENTER  XMIN  =  MIN  X  VALUE  AT  WHICH  TO  OUTPUT  CDF 

>  XMAX  =  MAX  X  VALUE  AT  WHICH  TO  OUTPUT  CDF 

>  DX  =  INTERVAL  OF  X  VALS  WHERE  CDF  IS  OUTPUT 

>  ALL  THREE  VALUES  ON  SAME  LINE;  E.G. 

> 

>  0.5  10.0  0.5 

> 

>  GIVES  OUTPUT  AT  20  X  VALUES  FROM  0.5  TO  10.0 

> 

>  ENTER  XMIN, XMAX, DX: 

5.0  20.0  0.5 

> 

>  ♦♦  ENTER  NFILES=NUMBER  OF  DATA  SET  FILES  TO  BE  FIT 

>  CURRENT  OPTIONS:  NFILES=  1  THROUGH  3 
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ENTER  NFILES: 


> 

> 

1 
> 

>  INPUT  FILE  NUMBER:  1 

>  ♦+  ENTER  FILENAME  WHERE  DATA  ARE  STORED, 

> 

>  ENTER  INPUT  FILENAME: 
gumbel . dat 

> 

>  ENTER  IDIST  =INDEX  OF  DISTRIBUTION  TYPE  TO  BE  FIT 

>  CURRENT  OPTIONS: 


> 

IDIST  = 

1  ... 

NORMAL 

> 

IDIST  = 

2  ... 

LOGNORMAL 

> 

IDIST  = 

3  ... 

EXPONENTIAL 

> 

IDIST  = 

4  ... 

WEIBULL 

> 

IDIST  = 

5  ... 

GUMBEL 

> 

IDIST  = 

6  ... 

SHIFTED  EXPONENTIAL 

> 

IDIST  = 

7  ... 

SHIFTED  WEIBULL 

> 

>  ENTER  IDIST: 

5 
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4  Output  Format  and  Extreme  Wave  Height 
Example 


The  first  output  section  of  FITS  provides  summary  statistics  for  each  of  the 
datafiles  considered.  These  include  (1)  sample  moments  from  the  data,  (2) 
predicted  moments  from  the  fitted  distribution,  and  (3)  underlying  distri¬ 
bution  parameters.  Comparison  of  (1)  and  (2)  can  serve  to  verify  the  fit  of 
low-order  moments,  and  the  agreement  of  higher  moments  (e.g.,  skewness 
and  kurtosis)  not  used  in  the  fitting  can  offer  a  rough  goodness-of-fit  test. 

The  second  section  of  output  gives  distribution  estimates  for  Gi{x),  the 
probability  that  a  future  outcome  exceeds  x  as  estimated  from  datafile  i.  It 
also  reports  the  distribution  of  annual  maxima,  Gann{.^')i  estimated  from 
Eq.  1.  Note  that  all  preceding  output  lines  begin  with  which  is  inter¬ 
preted  as  a  comment  within  the  public-domain  gnuplot  plotting  package. 
Thus  the  output  file  can  be  plotted  directly  with  gnuplot. 

The  contents  of  the  gumbell.out  output  file  follow  here.  The  output 
confirms  that,  as  intended,  the  fitted  Gumbel  model  preserves  the  mean 
m*=8.275m  and  standard  deviation  o-x=0.816m  estimated  from  the  data. 
The  Gumbel  model  tends  to  overestimate  the  higher  moments,  however;  its 
predicted  skewness  1.140  and  kurtosis  5.40  varies  notably  from  the  values 
-0.053  and  1.91  estimated  from  the  data. 

This  suggests  that  the  Gumbel  model  may  somewhat  overestimate  the 
chance  of  large  wave  heights  in  this  case.  This  has  also  been  observed  previ¬ 
ously  for  this  data  set  (Winterstein  and  Haver,  1991),  where  a  cubic  distortion 
of  the  Gumbel  model  was  introduced  to  capture  this  trend.  Here  we  seek 
to  model  this  effect  by  tail-fitting,  within  this  population  of  annual  maxiina. 
Specifically,  we  also  apply  FITS  to  the  same  data  with  a  shifted  exponential 
model  (IDIST=6).  Two  cases  are  considered,  corresponding  to  lower-bounds 
of  x/ot<;=8.0m  and  8.5m.  The  corresponding  input  and  output  files  in  these 
cases  are  gumbel?. in  and  gumbel? .out,  with  “?”=2  and  3.  The  input  files 
show  the  additional  input  line  required  for  shifted  distributions,  while  the 
output  (at  least  for  x/o^„=8.5m)  shows  a  somewhat  reduced  100-year  wave 
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{  Prob[  annual  max  Hs  >  h  ] } 


Figure  1:  Estimating  Annual  Max  Hs  distribution. 
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height  (10.4m  vs  10.8m  in  the  Gumbel  case).  These  distribution  fits  are 
shown  in  Figure  1.  They  suggest  that  these  shifted  models  can  begin  to 
capture  trends  in  distribution  tails,  not  available  in  moment-fits  to  the  entire 
population. 

A  listing  of  gumbel  1.  out  follows: 


# 

#  RESULTS  FOR  FILE  NUMBER:  1 

#  FILE  NAME  :  gumbel.dat 

#  DURATION  OF  DATABASE:  19.00 

#  DIST  TYPE  SELECTED:  GUMBEL 

#  NUMBER  OF  DATA  USED:  19 

# 

# 

#  MOMENTS  FROM  SAMPLE  DATA  (  MEAN,  SIGMA,  SKEWNESS,  KURTOSIS) 

#  0.8275E+01  0.8186E+00  -0.5282E-01  0.1905E+01 

# 

# 

#  MOMENTS  FROM  FITTED  DIST  (  MEAN,  SIGMA,  SKEWNESS,  KURTOSIS) 

#  0.8275E+01  0.8186E+00  0.1140E+01  0.5400E+01 

# 

# 

#  DISTRIBUTION  PARAMETERS  (SEE  DOCUMENTATION  FOR  DEFINITION) 

#  0.8275E+01  0.8186E+00  0.1567E+01  0.7906E+01 

# 

# 

#  **  FITTED  DISTRIBUTION  RESULTS  ♦♦ 

# 

#  PI  =  Prob  {Outcome  >  x}  in  database  1, 

#  P2  =  Prob  {Outcome  >  x}  in  database  2, 

#  .. . 

#  Pc  =  Prob  {Ann  max  >  x}  including  all  databases 

# 

#  X  PI  ... .  Pc 

0.5000E+01  O.lOOOE+01  0.6321E+00 
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0.5500E+01 
0.6000E+01 
0.6500E+01 
0.7000E+01 
0.7500E+01 
0.8000E+01 
0.8500E+01 
0.9000E+01 
0.9500E+01 
O.lOOOE+02 
0.1050E+02 
O.llOOE+02 
0.1150E+02 
0.1200E+02 
0.1250E+02 
0.1300E+02 
0.1350E+02 
0.1400E+02 
0.1450E+02 
0.1500E+02 
0 . 1550E+02 
0.1600E+02 
0.1650E+02 
0.1700E+02 
0.1750E+02 
0.1800E+02 
0.1850E+02 
0.1900E+02 
0.1950E+02 
0.2000E+02 


O.lOOOE+01 

O.lOOOE+01 

0.9999E+00 

0.9840E+00 

0.8489E+00 

0.5783E+00 

0.3260E+00 

0.1649E+00 

0.7905E-01 

0.3692E-01 

0.1704E-01 

0.7822E-02 

0.3581E-02 

0.1638E-02 

0.7486E-03 

0.3421E-03 

0.1563E-03 

0.7141E-04 

0.3263E-04 

0.1491E-04 

0.6810E-05 

0.3111E-05 

0.1421E-05 

0.6494E-06 

0.2967E-06 

0.1356E-06 

0.6193E-07 

0.2830E-07 

0.1293E-07 

0.5906E-08 


0.6321E+00 
0.6321E+00 
0.6321E+00 
0.6262E+00 
0.5721E+00 
0.4392E+00 
0.2782E+00 
0.1520E+00 
0.7600E-01 
0.3625E-01 
0.1690E-01 
0.7791E-02 
0.3575E-02 
0 . 1636E-02 
0.7483E-03 
0.3420E-03 
0 . 1563E-03 
0.7141E-04 
0.3263E-04 
0.1491E-04 
0.6810E-05 
0.3111E-05 
0.1421E-05 
0.6494E-06 
0.2967E-06 
0.1356E-06 
0.6193E-07 
0.2829E-07 
0.1293E-07 
0.5906E-08 
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A  Specific  Distribution  Types  and  Functional 
Forms 

Our  purpose  here  is  twofold.  First,  we  seek  to  define  the  precise  forms  of 
probability  distributions  available  in  FITS.  Second,  we  wish  to  indicate  pre¬ 
cisely  how  the  distribution  parameters  are  defined,  and  reported  in  FITS 
through  the  output  of  the  subroutine  DISPAR.  This  latter  information  may 
be  useful  if  the  user  seeks  to  perform  additional,  off-line  calculations  based 
on  the  fitted  distributions. 

The  following  basic  distribution  types  are  currently  available  within  FITS. 

IDIST=1:  Normal  Distribution.  The  cumulative  distribution  function  (CDF) 
is  given  by 

P[  Outcome  <  x]  =  F{x)  =  $( - -)  (2) 

O'x 

in  terms  of  the  standard  normal  CDF 

$(a:)  =  f  exp{-v,^ /2)du  (3) 

y/z'TT  -/-oo 

In  this  case  the  subroutine  DISPAR  returns  the  parameters  DPARM(l)=mj; 
and  DPARM(2)=a3;. 

IDIST=2:  Lognormal  Distribution,  The  CDF  is  given  by 

P[  Outcome  <  x]  =  F{x)  =  $( — — - ;  a:  >  0  (4) 

with  $(a;)  again  given  by  Eq.  3.  In  this  case  the  subroutine  DISPAR  re¬ 
turns  the  parameters  DPARM(l)=mx,  DPARM(2)=ax,  DPARM(3)=Tninx> 
and  DPARM(4)=(7ini. 

IDIST=3:  Exponential  Distribution.  The  CDF  is  given  by 

P[  Outcome  <  x]  =  F{x)  =  1  -  exp[-(-^)] ;  x  >  0  (5) 

TTtx 

In  this  case  the  subroutine  DISPAR  returns  the  parameters  DPARM(l)=mx 
and  DPARM(2)=ax,  and  no  additional  parameters. 
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IDIST=4:  Weibull  Distribution.  The  CDF  is  given  by 

P[  Outcome  <  x]  =  F{x)  =  1  -  exp[— (^)“] ;  x  >  0  (6) 

In  this  case  the  subroutine  DISPAR  returns  the  parameters  DPARM(l)=ma;, 
DPARM(2)=<Ta:,  DPARM(3)=q:,  and  DPARM(4)=^. 

IDIST=5:  Gumbel  Distribution.  The  CDF  is  given  by 

P[  Outcome  <  x]  =  F{x)  =  exp{—  exp[— a(x  —  u)]}  (7) 

In  this  case  the  subroutine  DISPAR  returns  the  parameters  DPARM(l)=ma;, 
DPARM(2)=(7a;,  DPARM(3)=a,  and  DPARM(4)=u. 

IDIST=6:  Shifted  Exponential  Distribution.  The  CDF  is  given  by 

P[  Outcome  <  x]  =  F(x)  =  1  -  exp[-(^ - ^)];  x>xiaw  (8) 

THx 

As  when  IDIST=3,  in  this  case  the  subroutine  DISPAR  returns  the  pa¬ 
rameters  DPARM(l)=ma:  and  DPARM(2)=cra:,  and  no  additional  pa¬ 
rameters. 

IDIST=7:  Shifted  Weibull  Distribution.  The  CDF  is  given  by 

P[  Outcome  <  x]  =  F{x)  =  1  -  exp[-( - ^  (9) 

As  when  IDIST=4,  in  this  case  the  subroutine  DISPAR  returns  the  pa¬ 
rameters  DPARM(l)=m:c,  DPARM(2)=a^,  DPARM(3)=a,  and  DPARM{4)=/?. 
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